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We show that a ferromagnetic Bose-Einstein condensate has not only hne-like vortex excitations, 
but in general also allows for point-like topological excitations, i.e., skyrmions. We discuss the thermo- 
dynamic stability and the dynamic properties of these skyrmions for both spin-1/2 and ferromagnetic 
spin-1 Bose gases. 
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I. INTRODUCTION 

An understanding of quantum magnetism is important for a large number of phenomena in physics. Three well- 
known examples are high-temperature superconductivity, quantum phase transitions and the quantum Hall effect. 
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f^ Moreover, it appears that magnetic properties are also very important in another area, namely Bose-Einstein conden- 

cn ■ 

[_^ ' sation in trapped atomic gases. This has come about because of two independent experimental developments. The 

f^ first development is the realization of an optical trap for atoms, whose operation no longer requires the gas to be 

C^ ' doubly spin-polarized ||J,p[- The second development is the creation of a two-component Bose-Einstein condensate 

S' . 

I ' PI , which by means of rf-fields can be manipulated so as to make the two components essentially equivalent B . As 

Ch ' 

Q ■ a result the behavior of both spin-1 and spin-1/2 Bose gases can now be experimentally explored in detail. Indeed, 

O ' 

►^ '■ at present already such diverse phenomena as domain walls ra] , macroscopic quantum tunneling p] , Rabi oscillations 
"— I . 

^% ' and vortices [Q have been observed. 

■ " " ' Theoretically, the ground-state structure of these gases has recently been worked out by a number of authors P-pd|] 



and also the first studies of the line-like vortex excitations have appeared P,n3Jl3l. However, an immediate question 
that comes to mind is whether the spin degrees of freedom allow for other topological excitations that do not have 
an analogy in the case of a single component or scalar Bose condensate. It is one of the aims of this paper to show 
that the answer to this question is in general affirmative. In particular, we show that a ferromagnetic Bose-Einstein 
condensate has so-called skyrmion excitations, which are nonsingular but nevertheless topologically nontrivial point- 
like spin textures. Roughly speaking, the skyrmion is an excitation that can be created out of the ground state, in 
which all the spins are aligned, by reversing the average spin in a finite region of space. Skyrmions are also known 
from nuclear physics iQ and the quantum Hall effect |p5|, but to observe them in an atomic gas would be exciting 



since in that case a completely microscopic understanding of their behavior is possible. In nuclear physics and the 
quantum Hall effect this is not true because of the nonpcrturbative nature of QCD and the presence of impurities 
that obstruct the center-of-mass motion of the skyrmions, respectively. Having proven their existence, we then turn 
to the investigation of the precise texture, the energetic stability, and finally the dynamical behavior of skyrmions. 
Some of the results of this analysis we have already reported in previous communications |16| , but here we try to give 
a much more complete and detailed picture of the skyrmion physics in spinor Bose-Einstein condensates. 

With this in mind, we would like to mention that we recently have also considered 't Hooft-Polyakov monopoles |17| ] 
in an antifcrromagnctic Bose-Einstein condensate |18|] . These topological excitations are in fact singular skyrmions, 
but have quite different properties than the nonsingular skyrmions which are the object of study in this paper. In 
particular, due to the singular nature of the spin texture of the 't Hooft-Polyakov monopole, the condensate density 
vanishes in the core and the monopole turns out to be thermodynamically stable. This is completely analogous to 
the case of a vortex in a scalar Bose-Einstein condensate. Both these features are not shared by the nonsingular 
skyrmions, which complicates the analysis considerably. The most important problem in this respect is that for a 
nonsingular skyrmion the topology allows for a spin texture with an arbitrary intrinsic size. As a result the stability 
of the skyrmion is now determined by energetic arguments and not by topological arguments as in the case of the 
singular 't Hooft-Polyakov monopole. 

The paper is organized as follows. In Sec. ||, we use the symmetry properties of the order parameter of a spinor 
Bose-Einstein condensate to show that quite generally skyrmion excitations indeed exist in such a condensate. We 
then turn our attention to the ferromagnetic case and discuss some general properties of these topological excitations, 



especially the spin texture, the superfluid velocity profile, and also the density profile. In Sec. [H we investigate 
the energetic stability of the skyrmion and show that from a thermodynamic point of view they are always unstable 
and tend to shrink to microscopic sizes. As a result we then consider the dynamical stability of the skyrmion and in 
particular determine the rate at which they collapse. We find that under certain conditions this rate is actually much 
smaller than the decay rate of the condensate itself due to various inelastic processes. Therefore, the skyrmion can 



for all practical purposes be considered as a (meta)stable excitation and we are justified to look in Sec. |V| also at 
other important dynamic properties of the skyrmion such as the "spin" and the center-of-mass motion. Finally, we 
end in Sec. |V|by a summary and some conclusions. 



II. SKYRMIONS AS TOPOLOGICAL EXCITATIONS 

In this section we discuss in detail the main static features of skyrniions. First of all, we show from the symmetry 
of the order-parameter space of the spinor condensate, that from a topological point of view point-like skyrmion 
excitations can indeed exist in both spin-1/2 and spin-1 Bose-Einstein condensates. Focussing then on the ferromag- 
netic case, we introduce a convenient parameterization of the skyrmion texture, which allows us to incorporate most 
easily the nontrivial winding number associated with the skyrmion. Next, we write down the energy functional for a 
ferromagnetic spinor condensate and, by substituting the above mentioned parametrization, derive the corresponding 
Euler-Lagrange equations for the density profile and the spin texture of the skyrmion. Finally, to simplify the actual 
calculation of the skyrmion density profile and the spin texture, we propose a variational approach that automatically 
takes into account the desired overall features of the skyrmion texture. 

A. Topological Considerations 

To find all possible topological excitations of a spinor condensate, we need to know the full symmetry of the 



macroscopic wave function ^(r) = ^/rl(r)({r), where n(r) is the total density of the gas, C(r) is a normalized spinor 
that determines the average local spin by means of (S)(r) = (* {r)S abCb{^), and S are the usual spin matrices obeying 
the commutation relations [6*0, Sp] = ieafj^Sj. Note that here, and in the following, summation over repeated indices 
is always implicitly implied. From the work of Ho [g] we know that in the case of spin-1 bosons we have to consider two 
possibilities, since the effective interaction between two spins can be either antiferromagnetic or ferromagnetic. In the 
antiferromagnetic case the ground-state energy is minimized for (S)(r) == 0, which implies that the parameter space 
for the spinor C(r) is only U{1) x S"^ because we are free to choose both its overall phase and the spin quantization 
axis. In the ferromagnetic case the energy is minimized for |(S)(r)| = 1 and the parameter space corresponds to the 
full rotation group 50(3). Using the same arguments, we find that for spin-1/2 bosons the order-parameter space of 
the ground state is always equivalent to SU{2) ||l9|l . 

What do these results tell us about the possible topological excitations [ popifl ? For line- like defects or vortices, we 
can assume C(r) to be independent of one direction and the spinor represents a mapping from a two-dimensional plane 
into the order-parameter space. If the vortex is singular this will be visible on the boundary of the two-dimensional 
plane and we need to investigate the properties of a continuous mapping from a circle S^ into the order-parameter 



space G, i.e., of the first homotopy group iti{G). Since ■ki{SU{2)) — T:i[SO{i)) = Z2 and 7ri([/(l) x S*^) = Z, we 
conclude that a spin-1/2 and a ferromagnetic spin-1 condensate can have only vortices with a winding number equal 
to 1, whereas an antiferromagnetic spin-1 condensate can have vortices with winding numbers that are an arbitrary 
integer. Physically, this means that by traversing the boundary of the plane, the spinor can wind around the order- 
parameter at most once or an arbitrary number of times, respectively. This conclusion is identical to the one obtained 
previously by Ho Q . 

If the vortex is nonsingular, however, the spinor C,{y) will be identical everywhere on the boundary of the two- 
dimensional plane and it effectively represents a mapping from the surface of a three-dimensional sphere S"^ into 
the order parameter space. We then need to consider the second homotopy group Ti2{G). For this we have that 
'K2{SU{2)) ^ 7r2(50(3)) — and tt2{U{\) x S*^) = Z. Hence nonsingular or coreless vortices are only possible for a 
spin-1 condensate with antiferromagnetic interactions. It therefore appears that the nonsingular spin texture discussed 
in Ref . ||§] , is topologically unstable and can be continuously deformed into the ground state by "local surgery" |gj . 

We are now in a position to discuss point-like defects. Since the boundary of a three-dimensional gas is also the 
surface of a three-dimensional sphere, singular point-like defects are also determined by the second homotopy group 
■K2{G). Such topological excitations thus only exist in the case of a spin-1 Bose gas with antiferromagnetic interactions. 
We call these excitations 't Hooft and Polyakov monopoles ||l^, although it would be justifiable to call them singular 
skyrmions. For nonsingular point-like defects the spinor Q[r) will again be identical on the boundary of the three- 
dimensional gas. As a result, the configurations space is compactified to the surface of a four-dimensional sphere S^ and 
we need to determine the third homotopy group 7r3(G). For this we find 'Kj,{SU{2)) ~ 'n^{SO{?>)) = tt3{U{1) x S'^) = Z. 
Hence nonsingular skyrmion excitations exist in all three cases. 

B. Skyrmion Texture 

We consider from now on only the case of a homogeneous and ferromagnetic spinor condensate. In the ground 
state all spins are aligned along the direction of a uniform and sufficiently weak magnetic field, which we take to be 
along the z-axis. The uniform magnetic field is needed only to direct the spins in the ground state, but it should not 
provide a substantial energy barrier for spin flips. The fact that we consider a homogeneous gas and not a confined 
one is only for simplicity and turns out not to be crucial for the practical applicability of our work. This is so because 



the, for our purposes relevant length scale over which the skyrmion spin deformations take place is always of the 
order of the correlation length, which under typical experimental conditions is much less than the length scale for 
density gradients due to the confining potential. The skyrmion excitation is a space-dependent spin deformation of 
the ground state and can thus be represented by a position-dependent spinor C(r)- A convenient way of introducing 
the position-dependence in the spinor is to write it in terms of a position-dependent rotation that acts on the constant 
spinor C^ associated with the ferromagnetic groundstate. In this manner we have 



C(r) = exp^--f^(r).SK^ (1) 



Here, the constant spinor C^ minimizes the Zeeman energy and is given by 

)' -^-2 



c^ = <: 



, 5"=! 

in the usual basis that diagonalizes S^- Furthermore, Jl(r) is a real vector function of r. It parametrizes the 
ferromagnetic order-parameter space, which due to our incorporation of the factor l/S* in Eq. (|l|) is always a sphere of 
radius tt. We point out that here and throughout the following S can only take the values 1/2 or 1. The significance 
of Eq. (|l|) is that at a point r the spinor C^ is rotated by an angle that equals |ri(r)|/S' around the direction of 
r2(r). There is no restriction on the generality of spin textures produced by this means; it is merely a convenient 
parameterization of the order-parameter space in terms of fi(r). Since we are mostly interested in the equilibrium 
properties, we assume here for simplicity the maximally symmetric shape of the skyrmion which is expected to have 
the smallest possible gradients. This means that we take 

r2(r) = LL>{r)vjr = w(r)f , (3) 

where the function ioir) should obey the boundary conditions (jj(0) — 27r and X\vi\j.^oa ^{f) — |^3l- Thus, at these 
boundaries the rotation operator in Eq. (^ becomes the identity and gives rise to C(r = 0) = C,(v — > oo) — C^ . 
Furthermore, w(r) should as a function of radius decrease monotonocally from 27r to 0, since this will correspond 
to the smallest gradient energy for the spin deformations. With this ansatz for r2(r) and its boundary conditions, 
we see that by traversing the whole configuration space, we exactly cover the order-parameter space twice, which is 
required to avoid a singular behavior of the spinor at r = 0. Indeed, the boundary condition w(0) — n, which in first 



instance appears to be the right one as it leads to a spin texture which covers the order-parameter space only once, 
is physically unacceptable because it results in a singular behaviour of the spinor C(r) and therefore in divergencies 
in the gradient energy. The reason for this singular behaviour is that, given the maximally symmetric shape of the 
skyrmion, the spinor C(r) is always equal to C^ on the z-axis. Note that apart from topological reasons, the condition 
linir^oo ^('") = also ensures that the spin deformations associated with the skyrmion have a finite range, so that 
only a finite amount of energy is required to excite the skyrmion from the ground state. Thus, both on the z-axis and 
far away from the center of the skyrmion the spins are directed as in the ground state. From these general features 
we can now draw a qualitative picture of the skyrmion texture. Suppose we approach the origin along a line through 
the origin, starting from a distance much larger than the range of the spin deformations. Then how does the average 
local spin behave as a function of distance? According to the above description of the rotation operator depicted in 
Fig. |lj, the spin initially points up along the z-axis, but by approaching the origin it will start to rotate around the 
radial direction until we approach the origin and the spin points up again. It is essential that the spins complete an 
integer number of cycles equal to Xj S in order to have a nonzero topological winding number given by 

which in this case is equal to 1. 

Summarizing, we have obtained the desired position dependence of C(r) in terms of a single function Lo(r) that 
describes how the average local spin vector is tilted from its orientation in the ground state. This function and the 
local density n(r) represent the two degrees of freedom that describe the skyrmion and our next task is to determine 
their precise spacial dependence. 

C. Energy Functional 

For a ferromagnetic spinor condensate. Ho ||8|] has shown that, within the mean-field approach, the energy functional 
in the absence of a magnetic field is given by 



E[n,C] = Jdr 



— (^VvMOj -/i7i(r) 



-1-— n(r)|VC(r)p + -T2^n2(r) 
2m 2 



(5) 



where m is the mass of the atoms, fi is their chemical potential and T^^ = Anah^ /m is the appropriate coupling 
constant that represents the strength of the interatomic interactions in terms of the positive scattering length a. The 
term ?i^n(r)|VC(r)P/2m represents the energy density associated with the gradients in the spin texture. It is this 
term that enforces the boundary condition lo{0) = 2tt, as discussed previously. 

For sufficiently large distances the gradient terms vanish and we infer from Eq. (^ that /i = UocT'^^ , where 
^oo — "(r -^ oo). Using this, Eq. (|g) can be put in a dimensionless form by scaling the lengths to the correlation 
length ^ = l/^/Svranoo, the density to Uoc, and the total energy to fPnao£^/2m. In this manner Eq. (|^) becomes 

e[fX] ^ j dp [(VpV7(p))' +/(p)|VpC(p)r + \f{p) f{p)\ , (6) 

with p = v/S, and n{v) ~ noof{p)- The gradient term | VpC(p)P is calculated in Appendix H explicitly as a function 
of uj{p). This is achieved by inserting in the usual basis the spin-1 and spin-1/2 matrices in the rotation operator, 
performing a power expansion and then using some properties of the powers of the spin matrices to sum the resulting 
infinite series into the following compact formula 

[(5-cos(20))(2HlM£)M) + 1(3 + COS (20)) (^) , S=l 
where is the azimuthal angle between p and the z-axis. 

Finally, the density profile /(p) and the function Lo{p) can be determined from the two coupled differential equations 

resulting from minimizing e[f, Q with respect to f{p) and uj{p), namely 



v7(p) 



|VpC(p)r + /(p)-l = 0, (8) 



and 



i<„„,.,(^)_..^ 



.'(/(.». (^)- 



0, (9) 



where 



i^Jdpfip), 5=1/2 

(/(P)>i - { (10) 

^ J dpi5- cos (26)), fip), 5=1 



and 



i^Jdpfip), 5 = 1/2 

(/(p)>2 = <; . (11) 

J-Jdpli3 + cosi20))fip), 5=1 
In the spin-1 case, Eqs. ([7|) and (|8|) show that the density profile is in principle anisotropic. However, the anisotropy 



turns out to be rather small as we show explicitly in Sec. |II 



D. Ansatz for uj{r) 

In principle, to calculate the exact skyrniion texture and density profile, we should now solve the two coupled 
nonlinear differential equations that result from minimizing the energy functional with respect to /(p) and C(p)j 
namely Eqs. (pf) and IM). As an alternative, we employ here a rather simpler, though less rigorous, approach. We 
simplify the calculations by introducing an ansatz for uj{p) that takes explicitly into account the physical boundary 
conditions discussed in the previous section. Our ansatz is 

u;ip)=4cot-'{p/X)\ (12) 

where A is a parameter that determines the radius at which Lo{p) crosses over from 27r to and physically corresponds 
to the size of the skyrmion. We see that this ansatz automatically satisfies the boundary conditions w(0) = 27r and 
linip^oo ^(/o) = 0, as required. It should be noted here that the detailed functional behavior of Lu{p) will turn out not 
to be crucial for our results on the skyrmion stability as long as it satisfies the prescribed boundary conditions and 
falls of monotonically. To substantiate this remark, we have performed also calculations using a number of different 
forms for uj{p) that satisfy the desired boundary conditions such as uj{p) = 2tt/{1 + [p/Xf'). We found, as expected, 
that the skyrmion energy only differs slightly from one function to the other and, in particular, that the energy is 
minimized for the function given in Eq. dl2). Furtunately, this ansatz is also simpler to handle analytically. 

III. STATIC AND DYNAMICAL STABILITY 

The most important question about the skyrmion excitation is whether it is energetically stable or not. In other 
words, how does the energy of the skyrmion depends on its size? As we show next, the answer to this question is 
that the skyrmion always tends to shrink to microscopic sizes to minimize its energy. Although this is an unfortunate 
result, it does not need to rule out an experimental observation of the skyrmion excitation as long as the typical 
time scale for this collapse is sufficiently long. Therefore, we also consider this problem after we have discussed the 
thermodynamical stability of the skyrmion. 

For sufficiently large skyrmions the gradients of the spinor Q{v) are small and density fluctuations are therefore also 
small. The energy of the skyrmion can then be approximated by {h'^noo/2m) J rfrj VC(r)p. If the size of the skyrmion 
is of order of A, this energy scales as A. This indicates that the skyrmion tends to shrink in order to minimize its energy. 



However, for smaller sizes the density fluctuations and their gradients start to grow and this simple argument no longer 
applies. For large skyrmions Eq. (g) shows that the density fluctuations scale as f{p) — f w |VpC(p)P oc f/A^ and 
their energy contribution thus behaves as f/A. Approaching smaller values of A, the kinetic energy term in Eq. (J^) 
increases and the density fluctuations will scale with a power of A that is greater than —2, because otherwise the 
density would become negative at some point. As a result the energy will scale with a power of A that is different 
from 1. Therefore, there is in principle a chance for stability if the energy associated with the density fluctuations is 
at some point increasing when A becomes smaller. To investigate this possibility we need to calculate the energy as a 
function of A exactly for all values of A. If the energy function possesses a (local) global minimum for a finite A, then 
the skyrmion is energetically (meta)stable. We have indeed calculated this energy curve, the details can be found in 



Sec. [[II B| below, and it turns out that the skyrmion energy actually increases monotonically with A. This means that 
a skyrmion of any finite size is energetically unstable and will tend to shrink to zero size. Of course, this condition 
holds within the Gross-Pitaevskii theory which describes only the long- wavelength physics. Corrections to this theory 
will lead to a finite, but microscopically small size of the skyrmion. However, as mentioned previously, for observing 
the skyrmion experimentally it is important to know the time scale for its collapse since it could be larger than the 
lifetime of the condensate itself due to various inelastic processes, such as two- and three-body collisions and, in an 
optical trap, photon absorption. 

From our analysis of the skyrmion shrinking rate, it turns out that there are two size regimes with different shrinking 
mechanisms. The first is for skyrmions with a size A that is much larger than the correlation length of the gas. We 
denote such large skyrmions as "adult" skyrmions. The second regime is for much smaller skyrmions with sizes of 
the order of the correlation length or less. We denote these small skyrmions as "baby" skyrmions. In the next two 
subsections we calculate the skyrmion energy and estimate its shrinking rate in these two regimes. 

A. Adult Skyrmions 

As mentioned earlier, adult skyrmions are large skyrmions with a size A that satisfies A 3> ^ = l/^/Siranoo, where ^ 
is the correlation length of a homogeneous gas of average density noc and s-wave interatomic scattering length a. In 
this section we express the lagrangian of the skyrmion in terms of a time-dependent skyrmion size A(i). An equation 
of motion for X{t) is then derived. It should be noted here that comparison between our calculation in this section 



and experiment is only meaningful for skyrmion sizes much less than the condensate size since this calculation is for 
a homogeneous system and docs not take into account the effect of the trap. The two conditions thus imply that the 
spinor condensate is deep in the Thomas-Fermi regime. 

For adult skyrmions the gradients in the spinor ({r) are small and therefore the fluctuations in the density Sn{r) = 
n{r) — rioo are also small compared to the average density riac,- The energy of the skyrmion can thus be determined 
by an harmonic approximation to the Gross-Pitaevskii energy functional given in Eq. (|5|). To second order in dn{r) 
it is given by 

^'["^C] = ^ Z^'" /*' Sn{r)x^\r - r')Sn{r') 

+ y'drn(r)(^|^|VC(r)p-7B.(S)(r)) , (13) 

where x(r ^ r') is the static density-density response function, which is defined by 

' -V^ + 167ranoo) xi^ - i"') = Hr - r') (14) 



Amric 
but explicitly reads 



muryo exp(-\/2|r-r ' | /£,) 
<ft2 |r-r'|/e 



^ir-r^) = ^ -^'\:Zur'"' ^ (15) 



B is either a fictitious, caused by resonant rf-fields, or a real homogeneous magnetic field, and 7 is the corresponding 
magnetic moment of the atoms in the trap. Considering first the ideal case of zero magnetic field B and solving for 
the density fluctuations induced by the spin texture, the lagrangian takes the form L[(\ = T[Q — E[Q, where 

EK]^ f dr n^^\VC{r,t)\' 
J 2m 

j dv j dv' |VC(r,t)|2x(r-r')|VC(r',t)|2 , (16) 



Sm? 



and 



drn^CHr,t)ih-ar,t) 



+ ljdvjdr' CHr,t)th^^C{r,t)x{r-r')CHr',m^^C{r',tl (17) 

as a result of the fact that in Gross-Pitaevskii theory, the action for the ferromagnetic condensate contains the time- 
derivative term J dt J dr n{r, t) (^ (r, t) ih d(^{r, t)/dt. The equation of motion for X{t) can be derived from the above 

10 



lagrangian by substituting our ansatz for ^(r, t) given in Eqs. (|i|), (H), and (|l^). In the present limit, X'3> ^ and the 
second term in the energy E[(^] is much smaller than the first one. Therefore, we do not take it into account. Using 
the explicit expressions for |VC(r)p fi-om Eq. (0), the energy is calculated to be 

E[X\ = c|-nooA, (18) 

zm 

where C = {29~20S)^/2n^. For the time-dependent contribution T[C] we calculate first C^(r, t)d(^{r, t)/dt which turns 
out to be equal to —iX{duj/dX) cos9 for both S — 1/2 and S — I, where the dot denotes a derivative with respect to 
time. Since duj/dX is an even function in r, the first term of T[Q vanishes. Thus r[A] turns out to be 

T[X] = -^mX\ (19) 

So combining Eqs. (|8|) and (|l^), we see that the action for the dynamical variable X{t) is equivalent to that of 
a particle with a position-dependent effective mass m* = \f2mXla in the linear potential V^(A) = Ch^nocX/2m. 
Scaling again A to the correlation length ^, the number y/2X/a takes the form X/^/4TTn^^a^. Experimentally, the 
dimensionless parameter rioofl^ is typically of order 10^'^. This leads to an effective mass to* « IO^Ato,. Thus a typical 
adult skyrmion mass, for say A = 10, is 10"^ times atomic mass. Finally, the equation of motion for X{t) reads 

2AA + A2 + c^0, (20) 

where c = (29 — 20S')7r. In this equation lengths are again scaled to the correlation length ^ and time is scaled to 
the correlation time r = 2mS,^/h. In these units the condition of validity of this equation is A 3> 1. For the initial 
conditions A(0) = Aq and A(0) = 0, we find that 



c 

^6 



A-Ao(l-^i2) (21) 



is a solution to Eq. (PQ). This formula shows that adult skyrmions shrink with a rate 7aduit ~ \/c/2Xot which 
indicates that spin-1/2 skyrmions shrink almost twice as fast as spin-1 skyrmions. In both cases one can make the 
shrinking rate smaller by exciting larger skyrmions. For realistic estimates of these shrinking rates, we restore the 
units in the above expression for 7aduit- For *^Rb spin-1/2 condensate of central density lO^^^cm^'^, the rate then 
reads 7aduit ~ 18. 06^/ A sec~^. 



11 



B. Baby Skyrmions 

The cruicial difference between adult skyrmions and baby skyrmions, apart from the size difference, is that density 
fluctuations for baby skyrmions are much more important than for the adult skyrmions. The density depletion 
produced by the spin gradients increases for smaller skyrmions and thus one can not consider the density to be 
essentially uniform anymore. Therefore, the linear-response approach followed in the previous subsection does not 
apply. To take the density fluctuations properly into account, we need to use the full energy functional. As already 
mentioned, it leads to a size dependence of the the skyrmion energy such that the skyrmion tends to shrink to zero 
size. The shrinking process, however, is now fundamentally different and more complicated than that for the adult 
skyrmions. In this case, while the skyrmion is shrinking the density becomes increasingly depleted in a region that 
forms a closed shell around the center of the skyrmion. For sufficiently small A the depletion will be so large that 
the atoms inside the closed shell are essentially isolated from the atoms outside the shell. This will considerably slow 
down the collapse of the skyrmion since the atoms within the shell need to tunnel over a potential barrier to escape to 
the other side of the shell and to enable the skyrmion to shrink further. It turns out that the size of the skyrmion at 
this stage is of the order of the correlation length and the lifetime of the skyrmion due to the above tunneling process 
can be much larger than the lifetime of the condensate itself. 

1. Density Profile and Texture 

With our ansatz for uo{p) in Eq. (n2[) the problem is significantly simplified since to determine the density profile 
we only have to solve Eq. (||), which is basically a nonlinear Schrodinger equation with some external potential 
|VpC(p)P- Substituting our ansatz for uj{p) in Eq. (0), the latter takes the form 

{onl P \2 [3+2(p/A)"+3(p/A)^] Q_i/o 

(22) 
4(17 + 3cos02)(^)2__i_^^ 5^1 

This is an off-centered potential barrier with a maximum height of 24.3/A^ at r sa 0.68A, for S — 1/2 and a maximum 
height equal to 3^3(17 + 3 cos^ 9) / AX"^ at r = A/Si/" for 5=1. Using this form for |VC(p)P we solve Eq. §) 
numerically for one particular value of A. Then we use the resulting density distribution f{p) to calculate the energy 
for that particular value of A. Performing the same calculation for different values of A we finally obtain the energy 
of the skyrmion as a function of A, from which we can judge the stability of the skyrmion. 

12 



For the spin- 1/2 case the density profile for various values of A is shown in Fig. || In the case of 5 = 1 the 
calculation is more complicated due to the angular dependence of |VpC(p)p. This complication is handled by 
noticing that |VpC(p)P can be rewritten in terms of Yoo{d,(f>) and F2o(^, 0) only, where Yim{0,<p) are the usual 
spherical harmonics. We thus also expand Eq. (||) in the spherical harmonics up to ? = 2 using 



V7(p) - yo{p)Yoo{e. c^) + y2ip)Y2o{e, </>). (23) 

Then two coupled equations can be obtained by taking the I = and I = 2 components of Eq. (^ . Specifically, we 
substitute Eq. ( p3|) into Eq. (g) and multiply by Foo(^i0) and Y2o{0,4>), respectively, and then perform an angular 
integration. The resulting equations take the form 

„2 16 (p/A)2 / 1 \ 1 3 3 2 n/5 3 

- ^> + A^ (i+(p/W r° ^ 7!"v ^ -y^ ^ ^'"'^ ^ ^y^ = '" ^''^ 

and 

72„ ^ 6 ^, ^ 16 (p/A)2 /I 58^^ \ 3 ^_ , ^ 3^/5^, ^ , ^ 15 . 



- ^,V. + -,y. + ^ (1 + (^/;,)4)2 (^7|2/o + y2/2J + ^2/2% + ^J/02/2 + ^y2 = V.. (25) 

where V^ = ^^ ( /^^^ ) • "^^"^ result of the numerical solutions of these two equations is presented in Fig. g. From 
this figure it is clear that the anisotropic part of /(p), which is represented by j/2(p), is considerably smaller than 
the isotropic part. This result will be employed in the following to simplify the energy calculation by neglecting the 
angular dependence of /(p). Note that for 5 = 1/2 this is not an approximation. 

Having specified Lo{r) enables also a detailed visualization of the skyrmion texture by calculating the average spin 
projections, i.e., < Sx >= C^(r)'S'2;C(r), < 5*^ >, and < Sz >■ In Fig. || we plot these quantities as a function of r for 
the spin-1 case using, again, A = ^. Similar plots for the spin-1/2 case can be found in our previous work |16[| . The 
skyrmion can be best visualized from its (5*2) (r) component, which is shown in Appendix |A| to be equal to 

(Sz) = cos^ 6* + cos w sin^ 9. (26) 

This expression is used to produce the two upper false-color figures in Fig. ^. These figures show clearly that the 
skyrmion corresponds to a torodial region. Along the inner radius of the torus, the direction of the average spin 
vectors is opposite to that in the ground state. An interesting property shows up when we calculate the associated 
superfluid velocity Vs(r) = — i?iC^(r)VC(r)/»Ti. This is given by 
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COS 6*—— f sinffsintj G sin6'siii^ i^/"^) 4> 

dr r r 



(27) 



From the above two expressions we find that the superfluid velocity is such that the atoms are simultaneously rotating 
around the inner radius of the torus and around the z-axis. This corresponds to a spiraling motion around the inner 
radius of the torus. The speed is maximum in the center of the torus. Furthermore, the maximum depletion of the 
density takes place at the maximum of |VC(r)P, which turns out to be in the center of the torus, where the spin 
vector is also completely flipped, i.e., (Sz) = —1. In the case of a spin-1/2 condensate, the spin texture consists of 
two torii. The velocity field in this case is such that atoms spiral around the center of one torus clock wise and along 
the center of the other torus counter clock wise, but for both torii the atoms rotate around the 2;-axis in the same 
direction. The spins are completely flipped along the inner radii of the two torii. The density is mostly depleted along 
a radius that is slightly smaller than the inner radius of the larger torus. The shift is only approximately 3 percent 
of the size of the skyrmion. 

2. Equilibrium State 

The outcome of the calculation in the previous subsection is the skyrmion density profile as a function of A. 
Therefore, the energy of the skyrmion, determined from Eq. (m, is now also a function of A. The skyrmion will be 
energetically stable if this energy E[X\ has a local or global minimum and the equilibrium size of the skyrmion will be 
equal to the value of A at that minimum. It should be noted that what wc call the energy of the skyrmion is actually the 
grand-canonical energy. This implies that the density profile is solved by using a fixed chemical potential for all values 
of A. Consequently, the number of atoms associated with the skyrmion excitation, which equals | J dr{n{r) — n.oo)|j 
is also A-dependent. (The subtracted term is to cancel the divergent part coming from the fact that the system is 
infinite. In a confined system there is no need to make this subtraction.) It may therefore not be immediately clear that 
minimizing the grand-canonical energy is appropriate. In principle, we must minimize the true energy of the skyrmion 
at a fixed number of atoms. We now explicitly show though, that these procedures are equivalent. Let us take the 
number of atoms in the system equal to N . This means that in the absence of the skyrmion J dr rioc = N. Suppose 
we now put in a skyrmion of size A in the condensate without affecting the asymptotic density, i.e., n(r — > cxd) = rioo- 
Due to the depletion of the density near the center of the skyrmion, the total number of atoms associated with this 
density profile is slightly less than N by an amount AiV(A) = j dr{noc — n{r,X)). Physically this means that to 
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produce this skyrmion we have to remove AiV(A) <C N atoms to the edge of the spinor condensate. If we, however, 
want to produce a skyrmion with the same number of atoms in the condensate, we have to adjust the asymptotic 
density. This impHes that a reasonable approximation to the density profile is actually {N/{N — AN))n{r, A) because 
then 

^ fdrnir,X)^N^^^4^^N. (28) 



N-AN J ' ' ' iV-AiV(A) 

Using this density profile to calculate the true energy of the skyrmion and expanding the various terms in it up to 
first order in AiV(A)/A^, we reproduce exactly the grand-canonical energy E[X]. Thus, we conclude that minimizing 
the canonical energy with a fixed number of atoms is indeed equivalent to minimizing the grand-canonical energy at 
a fixed chemical potential and thus with a changing number of atoms. The second method is clearly more convenient 
numerically since it is difficult to keep the number of atoms fixed for each value of A. 

In Fig. |5| we show our final result of this calculation, where we have taken n(r) to be isotropic, since in the 
previous section we saw that the angular-dependence of n(r) can indeed be neglected. In detail, we have plotted here 
E[n{r)] — iin{r) — {E[noo\ ~ M'^oo) against A. Again, we subtract the divergent part due to the infiniteness of the 
system. Unfortunately, both of the two curves do not have a minimum for any finite A and only contain a global 
minimum at A = 0. This means that, according to this figure, the skyrmion is energetically unstable, i.e., if the 
skyrmion is created with a finite size it will ultimately shrink to zero size. This process is clearly seen in Fig. g 
where we plot the skyrmion density profile for different values of A. In the next subsection we examine this shrinking 
process for the baby skyrmions more closely. In particular, we again look at the time scale over which this shrinking 
occures and find that for current experimental situations it is of the order of, or even longer than, the lifetime of the 
condensate itself. 

Finally, we mention here some remarks about the effect of a homogeneous external magnetic field and the slight 
differences in the scattering lengths of the different hyperfine spin states. In the presence of a homogeneous magnetic 
field B, it will be energetically unfavorable to flip spins. Therefore, this will lead to a further reduction in the spin of 
the skyrmion and, thus, not to a stabilization. A slight difference in the various scattering lengths Aa, as appropriate 
for the spin-1/2 ^^Rb condensate [£9|, leads to an additional term in the energy expression for the spinor condensate 
of the form {AnAati^ /m) J dr n{r){Sz){r), which has exactly the same effect as an external magnetic field. It can not 
stabilize the skyrmion either. 
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3. Shrinking Rate 

The goal of this section is to estimate the time scale over which the skyrmion shrinks to zero size. A key element 
in understanding the dynamics of this process is the behavior of the potential barrier produced by the spin gradients 
I VC(r)P when A goes to zero. This potential is written down explicitly in Eq. (p2[). When A is decreasing the location 
of the off-centered peak of this potential is approaching the origin as A and its height is increasing as 1/A^. This 
potential is thus an isotropic three-dimensional repulsive shell that is increasing in strength as it becomes smaller. We 
denote from now on the atoms inside this shell by the "core" atoms and the ones outside by the "external" atoms. We 
have seen in the previous subsection that, in principle, the skyrmion is unstable and that it is energetically favorable 
for the skyrmion to shrink to zero size irrespective of the initial size. Let us therefore consider a skyrmion of a large 
initial size such that the height of the corresponding potential barrier is less than the chemical potential /i = iiooT^^ 
of the system. The density distribution will be almost uniform everywhere except in the region where the barrier has 
its maximum height. In this region the density profile will be depleted only slightly. When the skyrmion starts to 
shrink the depletion of the density becomes larger, since the height of the barrier increases. This is achieved physically 
by transporting the core atoms over the potential barrier to the external region. So the central density n(0) decreases 
while the skyrmion is shrinking, as can be seen in Fig. 0. At a certain size the height of the barrier becomes equal 
to the chemical potential of the atoms at large distances. As the skyrmion shrinks even further the barrier height 
actually exceeds the chemical potential. At this stage the shrinking process slows down significantly since the core 
atoms will now have to tunnel through a potential barrier. The rate of this process will be characterized by the usual 
WKB tunneling rate, which we calculate now. 

To calculate this tunneling rate we assume that the skyrmion has shrunk to a size for which the barrier is so high 
that the overlap between the wave functions of the core atoms and the external atoms is exponentially small. The 
tunneling rate will then be determined by the chemical potential of the core atoms, as well as the barrier height. 
We consider the situation that there are Ncorc core atoms with a chemical potential /icorc less than the height of the 
barrier. Note that physically we must always have that /icorc ^ M since the collapse of the skyrmion will squeeze the 
atoms inside the core. For small r/\ we notice from Eq. ( p2[ ) that the potential (?i^/2m)| VC(r)P can be approximated 
by a harmonic potential well with a characteristic frequency ujq that equals \/48h/mX^ for S = 1/2 and \/32h/mX^ 
for S = 1. It thus has a characteristic length I — A/\/96 for S = 1/2 and I = A/8 for 5 = 1. This allows us to 
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use a Thomas-Fermi approximation, the vahdity of which will be discussed below, to calculate /icoro- Moreover, the 
tunneling rate is calculated using the following WKB expression 130] 



Wo 
7baby « TT- ^xp 
Ztt 



-2 f ' drV|VC|2-2m^corc/?i= 

J n 



(29) 



The radial points ri and r2 are the points where (fi.^/2m)|VC(r)p and ficorc intersect. The use of a Thomas-Fermi 
approximation is justified when the mean-field interaction energy of the core atoms is larger than the spacing between 
the lowest energy levels of the harmonic potential well. Specifically, the ratio 2Na/l = 2\/96Na/X should be bigger 
than 1. Figure |6|a shows the equilibrium size of the baby skyrmion, which is calculated by minimizing the total energy 
while keeping the number of the core atoms constant, as a function of the number of the core atoms A^coro- From 
this figure we observe that the above ratio equals approximately 1 for iVcorc = 4 and it increases for larger Ncoi-e- 
For experimental situations. Fig. |6|b shows that the skyrmion can live for a time ranging from seconds to hundreds 
of seconds, which is long enough to be observed. It should be noted that the above expression, which is based on 
a variational approach, somewhat overestimates the lifetime of the skyrmion. However, this is not crucial for our 
purposes, since by shrinking slightly the lifetime of the skyrmion increases considerably. 

IV. SKYRMION DYNAMICS 

A very important dynamical variable of the skyrmion arises from the fact that the Euler-Lagrange equations for 
the skyrmion spin texture is invariant under a space independent rotation of the average local spin (S) (r) around 
the magnetic field direction B |2J|. Mathematically, this means that if the spinor C'^'^lr) describes a skyrmion, then 
exp{— ii?B-S}C'''^(r) describes also a skyrmion with the same winding number and energy. The dynamics of the variable 
"(?(/:) associated with this symmetry is determined by the full action for the spin texture S[(^] ~ J dt{T[(^] — -E[C]), 



where E[(^] and T[(] are given by Eqs. (16) and (|17|) with Uoo replaced by the density profile n(r) of the skyrmion and 
x(r — r') should now be interpreted as the exact density-density correlation function in the presence of the skyrmion. 
Hence, substituting C{r,t) — exp{— ii?(t)B • S}C^''(r) and making use of the conservation of total particle number N 
to introduce the change of the average local spin projection on the magnetic field (AS'||)(r) = B- (S)(r) — iV^* induced 
by the skyrmion, we obtain, apart from an unimportant boundary term, that the dynamics of the rotation angle i!}{t) 
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is determined by the action 

where (AS'n"*) is the change of the total spin along the magnetic field direction and the "moment of inertia" of the 
skyrmion equals 

I = h'ldrldr' (A^||)(r)x(r,r')(A5||)(r') . (31) 

For an adult skyrmion this integral can be performed by substituting (A5||) = (Sz) — NS and using Eq. (jlj) for 
x(r,r')- The explicit expressions for (Sz) can be found in Appendix ^ An even simpler expression for x(r,r') can 
be optained by neglecting the gradient terms in Eq. (|lj). Then 

X(r,r') = ^<5(r-r'). (32) 

Using this approximation wc find, 



537r A"^ 
40V2 a '™' 


S^l/2 




S=l 



I^{ . (33) 
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The importance of this result is twofold. First, from the action in Eq. (30) we see that at the quantum level the 
hamiltonian for the dynamics of the wave function ^^(i?,^) becomes 



-^^(^-«<-r>)^ 



(34) 



Therefore, the ground state wave function is given by ^o(^) = e*^''/'\/27r, with K an integer that is as close as possible 
to (AS*,',"*). In this way we thus recover the fact that according to quantum mechanics the total number of spin-flips 
associated with the skyrmion texture must be an integer. More precisely, we have actually shown that the many-body 
wave function describing the skyrmion is an eigenstate of the operator B • S'°* with eigenvalue NS — K . Note that 
physically this is equivalent to the way in which "diffusion" of the overall phase of a Bose-Einstein condensate leads 
to the conservation of particle number |^,^ . 

Furthermore, the existence of this internal degree of freedom becomes especially important when we deal with more 
than one skyrmion in the condensate. In that case every skyrmion can have its own orientation and we expect the 
interaction between two skyrmions to have a Josephson-like contribution proportional to cos(^i — i?2)- As a result the 
phase diagram of a gas of skyrmions can become extremely rich and contain both quantum as well as classical, i.e., 
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nonzero temperature, phase transitions [|4|. In this context, it is interesting to mention two important differences with 
the situation in tlic quantum Hall effect. First, the fact that the spin projection K of the skyrmion is an integer shows 
that these excitations have an integer spin and are therefore bosons |27]] . In the quantum Hall case the skyrmions are 
fermions with half-integer spin, due to the presence of a topological term in the action 5[C] for the spin texture P8[ . 
Second, in a spinor Bosc-Einstcin condensate the skyrmions arc not pinned by disorder and arc in principle free to 
move. Both differences will clearly have important consequences for the many-body physics of a skyrmion gas. 

Focusing again on the single skyrmion dynamics, we can make the last remark more quantitative by using the ansatz 
C(r,t) = Cf^{Y — u(t)) for the texture of a moving skyrmion, which is expected to be accurate for small velocities 
du(t)/dt. Considering for illustrative purposes again only the isotropic approximation, we find in a similar way as 
before that the action for the center-of-mass motion of the skyrmion becomes 

^[u]^/rf4M(M^)\ (35) 

where the mass is now simply given by 

M=^JdrJdr' x(r, r') (v^) (r) • (v,) (r') (36) 

in terms of the superfluid velocity of the spinor condensate Vs(r). Note that in the anisotropic spin-1 case this mass 
is in principle a tensor, but with only very small nondiagonal components. 

Similar to the moment of inertia, we can calculate this mass explicitly for adult skyrmions. The result is 

M= ^ . (37) 

[ 2^A^, S=i 

The skyrmions thus indeed behave in this respect as ordinary particles. In contrast to Eq. (pG) there thus appears no 

term linear in du{t)/dt in the action ^[u]. This is a result of the fact that we have performed all our calculations at 

zero temperature. In the presence of a normal component, we anticipate the appearance of such a linear term with 

an imaginary coefficient. This will lead to damping of the center-of-mass motion of the skyrmion. It is interesting 

to note that if we perform the same analysis for a vortex in a scalar condensate, we find, due to the singular nature 

of the superfluid velocity field, even at zero temperature an additional term linear in the velocity of the vortex. In a 

quasi two-dimensional situation, it is in fact proportional to u(t) x du{t)/dt. This precisely results in the well-know 

Euler dynamics of vortices. 
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V. CONCLUSION 

One important conclusion of this work is that skyrmions in a ferromagnetic Bose-Einstein condensate are energet- 
ically unstable. However, we have also shown that the time scale on which the skyrmion shrinks can be of the order 
of, or even larger, than the life time of the condensate. It turns out that there are two very different mechanisms for 
the shrinking of the skyrmion. The first occurs for skyrmions with sizes much larger than the correlation length of the 
gas. In this case the shrinking rate is of the order of seconds for realistic experimental parameters. The second case 
concerns smaller skyrmions with sizes of the order or less than the correlation length of the system. In this case the 
shrinking rate is determined by the tunneling rate from the core of the skyrmion to the uniform part of the system. 
The tunneling takes place through a potential barrier developed by spin deformations. The typical time scale for such 
a shrinking process is of the order of 10 - 100 s. 

Although we have considered a homogeneous discussion, the inclusion of a trap in the system will not change the 
main conclusions of this work. The presence of a trap will have the effect that the uniform density rioo will be now 
a position-dependent quantity. Therefore, for our calculations to be valid for a trapped gas, the size of the skyrmion 
must be much less than the typical length scale for the density gradients of the confined condensate. Thus the 
condition A <C -R must be met, where R is the size of a spherical condensate. For typical experimental conditions 
we can use the Thomas-Fermi approximation to the size of the condensate |M|], which reads R = y^2/'K{Na/iy^^l 



where N is the number of atoms in the condensate and I — y^h/mujo is the characteristic length of a harmonic trap 
of frequency ajQ. Scaling R to the correlation length ^, the former takes the form R = 4:N^/^{l/a)^/^. For a *^Rb 
condensate with N — 10^ atoms and a trap frequency of the order of 27r x 100 rad/s this gives R w 500. This gives an 
estimate for the maximum size of a skyrmion such that the above condition remains valid. For example, our theory 
for adult skyrmions is applicable for the size range of A = 10^ up to say 50^. 

Another important remark about a confined gas is that skyrmions will tend to move towards the surface of the 
Bose-Einstein condensate since the density there is lower and thus their energy will be less. Apparently, this leads to 
a decrease in the shrinking rate since for both adult and baby skyrmions the shrinking time is proportional to 1/noo, 
see Eq. (pl|) and the discussion afterwards, and Eq. (29) and the discussion before it. However, our theory for the 



shrinking rates of skyrmions clearly breaks down on the surface since there the typical length scale for the density 
gradients due to the trap are of the same order as the correlation length. 
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Although skyrmion-antiskyrmion pairs can be created by the Kibble mechanism in a temperature quench or by 
sufficiently shaking up the condensate, a more controlled way of creating a skyrmion can be achieved by using a 
magnetic field configuration in which the fictitious magnetic field is always pointing radially outward and its magnitude 
increases monotonically from zero at the origin to a maximum value for large distances from the origin |32| . Applying 
this field configuration for such a long time that the spins at large distances have precessed exactly twice around the 
local magnetic field, creates a single skyrmion. Of course, for a real magnetic field the above configuration requires 
the use of magnetic monopoles, but for a fictitious magnetic field it can be achieved by appropriately tailoring the 
detuning, the polarization and the intensity of two pulsed Raman lasers. The required spatial dependence of the 
detuning can be created experimentally by separating the centers of the magnetic traps for the two spin species 
along the z-axis [Q. Furthermore, the desired behavior of the Rabi frequency can be achieved by making with the 
first Raman laser two standing waves in the x and y directions that are both polarized perpendicular to the z-axis. 
For the other Raman laser we only need to require that it produces a traveling wave with a polarization that has a 
nonzero projection on the z-axis, since we want to realize a Am — transition in this case. In the above geometry 
the skyrmion is created exactly in the plane where the detuning vanishes and in the nodes of the first Raman laser. 
Note that since the distance between these nodes is generally much bigger than the correlation length, we create in 
this manner a large skyrmion that will start to shrink but ultimately selfstabilizes at a smaller size. Once created, 
the skyrmion can be easily observed by the usual expansion experiments that have recently also been used to observe 
vortices |33| ]. Similar as with vortex rings, we then observe an almost complete depletion of the condensate in a ring 
around the position of the skyrmion. 
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FIGURE CAPTIONS 



FIG. 1. Schematic figure representing the action of the spin rotation operator for a maximaUy symmetric skyrmion. At a 
position r the average spin vector is pointing initially in the z-direction. The spin rotation operator rotates the spin vector 
around r by an angle u;{r). 



FIG. 2. The density profile for a skyrmion in a spin-1/2 condensate for different values of A. These curves are the numerical 
solutions of Eq. M) using Eq. (M) . 



FIG. 3. Density profile for a skyrmion in a spin-1 spinor condensate with a size of A = ^. The solid curve represents the 
isotropic part of n{r) and the dashed curve represents the anisotropic part. These are the numerical solution of Eqs. (124) and 



FIG. 4. False-color figures representing the average spin (S)(r) for a skyrmion in a spin-1 condensate in difi'erent Cartesian 
planes. The size of the skyrmion is taken here to be equal to A = ^. The (Sy) figures can be obtained from the (Sx) figures by 
using the axial symmetry of the texture. 



FIG. 5. Energy of a skyrmion as a function of its size. Plotted with the solid (dashed) curve is the energy for a spin-1 
(spin-1/2) spinor condensate. The dotted asymptotes are the large A limit of the energy, which is linear in A and is given by 
the spin gradient term in the energy functional only. Energy is in units of the correlation energy E^ = (47rft^/2?n^^). 



FIG. 6. Shown are (a) the baby skyrmions size as a function of the number of core atoms and (b) the shrinking rate of 
skyrmion as a function of the number of core atoms. This calculation was performed for a *^Rb spin-1/2 condensate with a 
scattering length of a = 5.4 nm. 
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APPENDIX A: DETAILS OF THE CALCULATION OF {So:), {Sy), {S,), C'''(r)VC(r), AND IVCp(r) 



We start by expanding the exponential operator in Eq. (nh, taking into account the simphfication expressed in Eq. 
(0). Inserting the exphcit forms of the spin-1/2 and spin-1 matrices given by [p3 



1 

M 1 



-1 
2 M 



1 
-1 



z, 



s = <^ 



010\ /0-10\ /lOO 

^(101 U+^ 1 -1 y+ 
^U 1 0/ ^l 1 / I -1 



s-h 



S^l 



(Al) 



and noting that (S • f )'^ equals S • r for A; =odd, and equals (S • f )^ for k =even enables a resuming of the odd and even 
powers of the exponent separately. As a result, one can put this exponential operator in the following simple form 



(A2) 



exp <; ~^^{r) ■ S| = exp |-^c^(r)f • S 

lcos(w(r))-i(f •S)sin(tj(r)), S=^ 

1 - (f • Sf [I - cos (cj(r))] - i{r ■ S) sin (w(r)), S = 1 
where 1 is the identity matrix. It is then straightforward to derive the expressions for the important physical 
quantities such as the average spin components (Sx) = C^{r)SxC{''^), (Sy) and (Sz), and the superfluid velocity 
Vs = — i^C^(r)VC(r)/?Ti. For S = 1/2 and 1, these quantities take the form 



(S*^) /5' = sin 6* sin sin {uj/S)+ sin (261) cos {(/>) sin^ (w/25') , 



(A3) 



{Sy)/S = - sin 6* cos ^ sin (w/5) + sin {29) sin (0) sin^ (w/25') , 



(A4) 



{Sz)/S = cos2(6i) + cos (u/S) sin^ 



(A5) 



and 



= ft 



7 Q O C* 

cos 6*-- f sin 6* sin {uj/S) 6 sin ^ sin^ (tj/2S') 4> 

dr r r 



We can now derive the expression for the spin gradient energy h'^\'VC{r)\'^ /2m as Eq. ( p2| ) shows apart from the 
constant fi^ /2m. 
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This figure "figure l.j peg" is available in "jpeg" format from: 



http://arXiv.org/ps/cond-mat/0107302vl 
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